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Precision measurements of SM Higgs recoil mass and cross-section are performed with 
the Higgs- Strahlung process e + e~ — > Z°h°, with Z° — » n + fj,~ and Z° — » e + e~ at 230 
GeV and 250 GeV center of mass energies, presuming Mh = 120GeV, assuming the 
integrated luminosity is 500/6" , based on full simulations of several detector models 
for the ILC. Results are given and compared in this article, of the two center of mass 
energies, and of several detector models. According to the studies, the Higgs mass 
can be determined with a statistical error of 39 MeV and 3.6% statistical error for the 
cross-section. 

1 Introduction 

The Higgs-strahlung process e + e~ — > Z°h° offers an unique opportunity for a model in- 
dependent precision measurement of the Higgs Boson mass by means of the recoil mass to 
the Z°, gives as M? = s + M| — 2Ez"y/s, where Z° decays to e + e~ (e channel) or /i + \i~ 
(fj, channel). At the same time, the Higgs production cross-section and therefore also the 
coupling strength at the Z°h° vertex can be determined, g 2 oc a = N/Ce. 

This article presents results of the studies performed for two center of mass energies 
(\/s), at 230 GeV and 250 GeV. The 230 GeV is the optimal y/s in terms of the resolution of 
Higgs Mass according to a previous study [H], where the e channel is analyzed using detector 
model LDC01Sc[2] (LDC1). The 250 GeV is the benchmark scenario as proposed in the 
ILD[4] LOI, where both the e channel and \i channel are analyzed using detector models 
LDCPrime_02Sc[2] (LDCP), LDC01_06Sc[l] (LDC6) and LDC_GLD_01Sc[2] (LDCG). 

2 Experimental Remarks 

In this study, the Higgs mass is presumed to be 120 GeV, the luminosity is assumed to be 
500 fb , and unpolarized beams are assumed. The events generator is PYTHIA[5], simula- 
tion is performed using MOKKAJ5], and the reconstruction is performed using MarlinReco[Sj 
and PandoraPFA[7]. Both of the two yfs analyses include Bcamstrahlung, ISR and FSR. The 
Beamstrahlung spectrum is generated using GUINEA-PIG [8], and the interface to PYTHIA 
is CALYPSO 0. 

The cross-sections for signals and backgrounds are shown in Tab. [T] For \fs at 230GeV, 
since only detector acceptance cut is applied, the Bhabha scattering (e + e~) gives a cross- 
section five orders of magnitude larger than that of the signal. The four detector models 
under study are composed of: TPC, time projection chamber; VXD, 5 single or 3 double 
layers of vertex detector; SIT/SET, 2 cylindrical layers of silicon strips inside and outside the 
TPC; FTD, pixel silicon disks in the forward region; ECAL, SiW electromagnetic calorime- 
ter; HCAL, scintillator hadronic calorimeter. The magnetic fields in tracking system are 4 
Tesla for LDC1 and LDC6, 3.5 Tesla for LDCP, and 3 Tesla for LDCG. The momentum 
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yfs = 2'AOGeV e channel 


Reactions 


er 


ZU^u ^ e + e -x 


6.3 fb 




5.96 x 10 b fb 


t + t~ — > e + e~Av 


146 fb 


W+W~ -> e + e~2v 


181 fb 


Z u h*Z v h* -> e+e-ff 


113 fb 



yfs = 250GeV e and fj, channels 


Reactions 


er 


Z"h u -> e+e-X 


7.5 fb 


Z u /i u -> fi+ii-X 


7.5 fb 


Z U Z U -> e+e"// 


78.7 fb 


Z U Z U -» 


79.0 fb 



Tabic 1: Signals (bold) and backgrounds, and their cross-sections. Note, for \fs 230GcV 
backgrounds, only detector acceptance pre-cut is applied, \cos8 e ± | < 0.98. Assume C = 
500 fb and unpolarized beams. 



resolution (AP/P 2 ) is about 5 x 10 5 GcV 1 for the four detector models for momentum 
within 10 to 100 GeVQU]. 

A cut based electron ID is developed for the analysis at y/s = 230 GeV. The efficiency 
is greater than 99.5%, and the rejection rate is 100% for /i, and larger than 98% for n, for 
momentum within 20 to 80GeV, which covers the momentum range of the e + /e~ candidates 
of the signal. For the y/s at 250GeV analyses, MC information is employed temporarily for 
the lcpton ID. 

The final state lepton pair selection is based on a \ 2 criteria that the invariant mass is 
the closest to the Z° mass. 



3 Background Rejection 

Standard Model (SM) Higgs with M h ~ 120 
GeV dominantly decays to 2 jetspj]. To- 
gether with e + e~ decayed from Z°, more 
than 95% of the signal have a multiplicity 
of larger than 4 in the final states. For the 
backgrounds Bhabha, tt and WW, in their 
visible final states, the multiplicity is 2. Due 
to bremsstrahlung, one electron may lead to 
several reconstructed tracks. Therefore, for 
the y/s at 230 GeV e channel analysis, cut 
on Ntracks > 6 is applied to reject all the 
three backgrounds totally, with a penalty of 
~ 8% reduction of the signal. The N tr acks 

distributions of signal and the three backgrounds are shown on Fig [T] 

A Likelihood method is applied for the rejection of ZZ background for both center of 
mass energies. 

The Likelihood of an event to be the signal is defined as Ls = Y\Pf, where the Pf is 
the probability of the event to be the signal according to the PDF of the signal of the ith 
selection variable. Similarly, the Likelihood of an event to be the background is defined as 
Lb = Y[PP ■ Thereafter, the Likelihood Fraction is defined as Jl = Ls/(Ls + Lb), which 
is within (0, 1). 



Figure 1: N tra cks distributions of: ZH top- 
left; Bhabha top-right; tt bottom-left; WW 
bottom- right. 
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Figure 3: Number of signal and number of back- 
grounds vs. Jl cut. 
cuts (bottom). 



(top), and S/y/WTB vs. f L 



In order to be less dependent on 
the physics presumptions and to ob- 
tain a flat background distribution af- 
ter its rejection, in this study, only two 
angular variables are employed: acol- 
linearity and cos9 e - vs. cos9 e + . The 
PDFs of these two variables of yfs 230 
GeV e channel are shown in Fig. O 

The number of signal and back- 
ground vs. Jl cuts, and the signif- 
icance S/\/S + B vs. Jl cuts are 
shown in Fig. [3j which have the N tra cks > 6 cut applied already, and they are determined 
within the fitting range of M/j from 118 to 135 GeV. 

According to the maximum of the sig- 
nificance, the fi cut is determined to be 
/l > 6.7, resulting in an overall signal se- 
lection efficiency of 55.5% and a fraction of 
remained ZZ background of 3.5% for y/s at 
230 GeV analysis. For the y/s at 250 GeV 
analysis, since only ZZ background is con- 
sidered, the same Likelihood method is ap- 
plied for the rejection of ZZ as that of the 
yfs at 230 GeV. The resulted efficiencies are shown in Tab. [2j 



Detectors 


LDCP 


LDC6 


LDCG 


e+e-X 


39.5% 


37.2% 


39.5% 




55.5% 


52.7% 


35.5% 



Table 2: Efficiencies of signal selection of yfs 
250 GeV e and \i channels, determined within 
M h window 118 - 135 GeV. 



4 Fitting Methods 

Referred to fitting methods in previous studies [12] [13] [14] |15j [16j . the Gaussian core for the 
Peak with Exponential complement for the Tail (GPET) method is chosen and improved in 
this study, which eventually has both it-self and its first derivative continuous at all points, 
as shown in EqJT] It is a partial function, the left part is a pure Gaussian, the right part 
is a sum of Gaussian and Exponential with the fractions of contribution to be (3 and 1 — f3 
respectively, and a factor k is introduced in order to keep the maximum (xo) been covered 
by the pure Gaussian. Fig[4] showed an example fitting to the signal using this formula with 
a x 2 /Ndf ~ 1. 
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f(x) = N 



(3e 



x^xo < f, 
a — 

> k 



(1) 



Signal Filling: 1631 H 



bela = 0.46 ±0.030 
k= 0.057 ±0.0030 
meanO = 1 20.000 ± 0.01 4 
Shifts 0.25 ±0.023 
slgma - 0.33 ±0.015 




Figure 4: Example fitting using 
the GCET formula. 



The GPET fit provides an explicit description of the 
final spectrum including the detector response. Together 
with the Higgs mass, the mass resolution can also be 
measured by this method. 

However, due to the uncertainty of the detector re- 
sponse, the maximum of final recoil mass spectrum has 
a shift to value larger than the physics presumption 
(M^ c = 120 GeV in this study). A correction of the de- 
tector effects is required in order to restore to real mass. 
Since the shift (A.t = xq — M™ c ) comes from the detec- 
tor effects, an explicit measurement of this shift from full 
detector simulation is possible and acceptable. 

For ^/s at 230GeV analysis, beside M™ c = 120 GeV, 
simulations and reconstructions for M™ c = 117, 118, 

119, 121, 122 and 123 GeV are performed, too. The shifts (Ax) are determined subse- 
quently. The mean of the shifts measured Ax = Q.27i)GeV is taken as the correction, and 
the standard deviation <r(Ax) = 0.021GeV, which is the uncertainty of the measurement 
of the shifts, is taken as the systematic error, for the Higgs mass measurement. For s/s 
250GeV analyses, since no simulations are available for M™ c besides 120GeV, the shifts 
measured from M™ c = 120GeV simulations are taken as the corrections, while the statisti- 
cal errors from the measurements of the shifts are taken as the systematic errors. Actually, 
this systematic error can be eliminated by increasing the statistics of MC data samples in 
the determination of the shifts. 

To describe the backgrounds spec- 
trum, the Chebyshev Polynomial 
function with two coefficients is em- 
ployed. 



Higgs Recoil Mass Spectrum, E =230GeV, LDC01SC 



5 Results 

The final fitting to the signal with 
background is performed with the 
mass (M/i = x — Ax), the mass res- 
olution (er in the pure Gaussian) and 
the number of signals as free parame- 
ters. A typical result is shown in Fig{5j 
and the fitting results including the 
systematic errors of Mh are listed in 
Tabll 



• Data 
— Signal 

Backgrounds 




Figure 5: The recoil mass spectrum of signal and 
backgrounds of y/s 230GeT^ e channel. 
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(GeV) 


Detector 
Model 


Chan- 
nel 


M h (GeF) 
(±siai. err. ± sys. err.) 


o ifb) 
(±stat. err.) 


S m (MeV) 
(±stat. err.) 


230 


LDC1 


e 


120.022 ± 0.039 ±0.021 


6.41 ± 0.43(6.7%) 


360 ± 17 


250 


LDCP 


e 


119.973 ±0.047 ±0.039 


7.82 ± 0.52(6.6%) 


540 ± 25 




120.019 ± 0.023 ±0.016 


7.78 ± 0.28(3.6%) 


500 ± 12 


LDC6 


e 


119.963 ±0.047 ±0.044 


7.93 ± 0.49(6.2%) 


560 ± 28 


/" 


119.994 ± 0.023 ±0.016 


7.45 ± 0.27(3.6%) 


550 ± 12 


LDCG 


e 


119.973 ±0.051 ±0.044 


7.24 ± 0.52(7.2%) 


490 ± 27 


A* 


120.003 ± 0.029 ±0.020 


7.45 ± 0.32(4.3%) 


530 ± 15 



Table 3: Results of Higgs mass (M h ), cross-section (a) and mass resolution (S m ). 



6 Conclusion and Outlook 

According to the analysis results shown in TabJ3J with the same luminosity, the statistical 
error of /x channel is about one half compared with the e channel. This is due to the 
bremsstrahlung effect of electrons reduces the statistics on the maximum of the recoil mass. 
The results between different detector models are nearly the same, since the momentum 
resolutions are roughly the same[T0]. 

Improvements of this analysis can be expected with polarized beams, which may increase 
the cross-section of Higgs-strahlung reaction, and the left-handed polarized positron beam 
with right-handed polarized electron beam may suppress the WW background largely. For 
the rejection of Bhabha scattering and ee — > fifx in a model independent analysis, the Pr of 
ISR photon can be used to balance the Pr of the di-lepton system. Taken the Higgs mass to 
be 120 GeV, the ZZ background can be further reduced by reconstruction both Z bosons. 
Finally, the fitting range and fitting method may need to be optimized. 

References 

[1] Presentation: 

http: //ileagenda. linear collider . org/contributionDisplay . py?contribld=405&sessionld=16&conf Id=2628 

[2] P. M. Freitas, MOKKA http://mokka.in2p3.fr/ 

[3] F. Richard, P. Bambade, ELAN Document-2007-001 (2007) 

[4] ILD (International Large Detector concept) Ihttp:/ /www.ilcild.org/l 

[5] T. Sjostrand, S. Mrenna, P. Skands, PYTHIA 6.4, JHEP 05 (2006) 026 

[6] ILC Software http: //ilesoft .desy.de/ 

[7] M. Thomson, AIP Conf. Proc. 896 (2007) 215-224 

[8] D. Schulte, PhD Thesis, University of Hamburg, (1996) 

[9] D. Schulte, CALYPSO: A Library to Include the Luminosity Spectrum in Event Generators (2006) 

[10] S. Aplin, A. Raspereza, LCWS08, Tracking Studies for ILD, (2008) 

[11] A. Djouadi, J. Kalinowski and M. Spira, Comput. Phys. Commun. 108, 56 (1998). 

[12] J-C. Brient, LC-PHSM-2000-049 (2000) 

[13] P. Garcia-Abia, W. Lohmann ,A. Raspereza, LC-PHSM-2000-062 (2000) 

[14] P. Garcia-Abia, W. Lohmann, Eur. Phys. J. direct C2:2, (2000) 

[15] M. Ohlerich, W. Lohmann, A. Raspereza, A. Schalicke, LCWS07 Proc. Vol. 1, 159-165 (2007) 

[16] M. Ruan, Y. Gao, F. Richard, Z. Zhang, ELAN Document-2008-010 (2008) 



LCWS/ILC2008 



5 



